Online Correction of Dispersion Error in 2D Waveguide Meshes 
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Abstract 

An elastic ideal 2D propagation medium, i.e., a mem- 
brane, can be simulated by models discretizing the 
wave equation on the time-space grid (finite-difference 
methods), or locally discretizing the solution of the 
wave equation (waveguide meshes). The two ap- 
proaches provide equivalent computational structures, 
and introduce numerical dispersion that induces a mis- 
alignment of the modes from their theoretical posi- 
tions. Prior literature shows that dispersion can be 
arbitrarily reduced by oversizing and oversampling the 
mesh, or by adopting offline warping techniques. In 
this paper we propose to reduce numerical dispersion 
by embedding warping elements, i.e., properly tuned 
allpass filters, in the structure. The resulting model ex- 
hibits a significant reduction in dispersion, and requires 
less computational resources than a regular mesh struc- 
ture having comparable accuracy. 

1 INTRODUCTION 

Membranes are the crucial component of most percus- 
sion instruments. Their response to an excitation, and 
their interaction with the rest of the musical instru- 
ment and with the environment, strongly affect the 
sound quality of a percussion. Physical modeling of 
membranes has drawn the attention of the computer 
music community when a new model based on the 
Digital Waveguide was designed, called 2-D Digital 
Waveguide Mesh |Van Duync and Smith, 1993fl . The 
model was proved to provide a computational struc- 
ture equivalent to a Finite Difference Scherme (FDS). 
In particular, it was shown that the numerical arti- 
facts introduced by the model cause a phenomenon 
called dispersion. This means that, even in a 
flexible medium, different spatial frequency compo- 
nents travel at different speeds, a nd this speed is 
direction- and frequency-dep endent [ ritrikwerda, 1989 , 
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Different mesh geometries have been stud- 
ied: each of them have an equivalent FDS, 
and exhibits its peculiar dispersion error func- 
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dispersion error is, with good approximation, in- 
dependent from the direction of propagation of 
the spatial components; at the same time, the 
Triangular Waveguide Mesh (TWM) defines, from 
a signal-theoretic point of view, the most effi- 
cient sampling scheme among the geometries that 
can be derived from mesh models used in prac- 
tice 
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Fontana and Rocchesso, 2000 1 . 
The independency from direction has been success- 
fully exploited pavioja, 2000 1 to warp the signals 
produced by the model, using offline filtering tech- 
niques [ Harma et al., 2000| . In this paper we work 
on a similar idea, but warping is performed online 
by cascading each unit delay in the TWM with a 
first-order allpass filter. By properly tuning the filter 
parameter, we will prove that a considerable reduction 
of the dispersion error can be achieved in a wide range 
around dc. 

This result is then compared with the performance of 
a TWM, oversized in order to reduce dispersion in the 
first modes. It will be shown that the warped mesh is 
less expensive in terms of memory and computational 
requirements. This evidence holds both for the straight 
waveguide and the FDS implementations. Our conclu- 
sion is that the most efficient, low-dispersion compu- 
tational scheme for membrane modeling is a triangular 
FDS where the unit delays are cascaded with properly 
tuned allpass filters. 

Having an efficient and accurate membrane model 
is a key step toward the construction of affordable, 
tunable, and realistic models of complete percussion 
instruments. In particular, the coupling bet ween 
air and membrane [ Fontana and Rocchesso, 1998 1 , and 
the interface with resonating structures are funda- 
mental components that should be added to the 
membrane model in order to achieve better real- 
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tion [Fontana and Rocchesso, 2000 



geometry exhibits two valuable properties: the 



2 ONLINE WARPING 

For a wave traveling along the waveguide mesh, the 
numerical dispersion error is a function of the two 
spatial frequency components. In the TWM, this 
function is symmetric around the origin of the spa- 
tial frequency axes, with good approximation. Con- 
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Figure 1: Plot of the dispersion error versus temporal 
frequency magnitude 



sequently, it makes sense to plot the dispersion fac- 
tor as a single-variable function of spatial frequency, 
moving from the center of the surface to the ab- 
solute band edge along one of the three directions 
defined by the waveguide orientations^ Assuming 
the waveguides to have unit length, the spatial band 
edge results to be equal to 2n/^/3 [rad/spatial sam- 
ple] [ Fontana and Rocchcsso, 2000| . A plot of the dis- 
persion factor versus temporal frequency is then cal- 
culated recalling the nominal propagation speed fac- 
tor, equal to l/v2 [spatial sample / time sample], af - 
fecting any finite difference model |5trikwerda, 1989 |. 
that fixes the edge of the temporal frequency at the 
value \/2ir/y/3 [rad/sample]. Figure |] depicts a plot of 
the dispersion factor. 

This analysis is confirmed by simulations conducted 
over a TWM modeling a square membrane of size 
24 x 24 waveguide sections, clamped at the four edges, 
excited at the central junction by an impulse. In 
fact, the impulse response taken at the central junc- 
tion shows that the positions of its modes match well 
with the theoretical frequencies of the odd modes of 
the membrane, each one of them being shifted by 
its own dispersion and by the nominal propagation 
speed factor. The results are depicted in Figure g, 
where the frequency response of the model is plot- 
ted together with (o) the theoretical positions (com- 
pressed by the nominal propagation speed factor) of 
the modes resonating in the membrane below \/2ir/y/3 
[rad/sample], and (x) the real positions of the same 
modes, affected by dispersion. Overall, dispersion in- 
troduces a modal compression that increases with fre- 
quency. The careful reader will note a slight difference 
between the calculated frequency cut and the band- 
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Figure 2: Frequency response taken at the center of 
a TWM (size 24 x 24) excited by an impulse at the 
same point. Theoretical positions of the odd modes 
resonating in a mebrane below v / 2""/v / 3 [rad/sample], 
weighted by the nominal propagation speed factor (o). 
Positions of the same modes affected by dispersion ( x ) . 



width of the signal coming from the simulation. This 
difference is probably due to the simplifying assump- 
tion of considering the dispersion function as direc- 
tion independent. Moreover, some modes show up as 
twin peaks. This may be due to the actual irregular 
shape of the resonator model, caused by the impos- 
sibility to design a perfectly square geometry using a 
TWM model. In order to conduct a controlled ana- 
lytical study we avoided using interpolation along the 
edge, even though this is recommended in practical 
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implemcntations [Aird et al., 2000 1. 
Let H (z) be the transfer function of a TWM, regard- 
less of the excitation (input) and acquisition (out- 
put) positions. The transfer function can be han- 
dled by conformal mapping, a method consisting in 
the application of a particular map T to the z- 
domain, in order to obtain a new, warped domain 
I = T(z) jMoorer, 1983| , |Harma et al., 2000| . The fre- 
quency response H(e JUJ ), calculated from H(z)\ s=e jQ , 
changes according with the map. 

Practical implementations of transfer functions ob- 
tained by conformal mapping are often affected by 
non computable loops, that can sometimes be re- 
solved Harma, 1998 1. In a TWM, delay free loops ap- 
pear whenever the map does not allow to extract an 
explicit unit delay. However, if we change the number 
of unit delays in each waveguide section of a waveg- 
uide mesh, we only change the number of Fourier im- 
ages of the frequency response^, by simply compressing 
each single image. Now, imagine a map that translates 
each unit delay into the cascade of a first-order allpass 



tude around the origin is constructed, resulting in a slight dif- 
ference respect to the curve adopted here. 
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time linear filter. 



is applied to a discrete- 
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Figure 3: Mapping functions z^ 1 — z~ x A{z) for 
equally-spaced values of the parameter a of the allpass 
filter A(z). Top line: a = 0. Bottom line: a = —0.9. 



filter A(z) and a unit delay, z^ 1 = z^ 1 A(z). If the 
allpass filter has a negative coefficient, the phase de- 
lay introduced by the filter ranges from one sample (in 
high frequency) to a certain value larger than one (in 
low frequency). Therefore, it is reasonable to expect 
an extra image (due to doubling the unit elements for 
high frequencies) , and a degree of compression that de- 
creases with increasing frequency. This is exactly the 
kind of behavior that is desired in order to counterbal- 
ance the effects of numerical dispersion. In order to 
get rid of extra frequency components it is sufficient 
to lowpass filter at the desired cutoff frequency, and to 
restore the correct positions of low-frequency partials 
it is sufficient to increase the temporal sampling rate. 
The frequency domain is warped according with the 
formula 



Q = arctan 



2 sinw{a + cosw} 



a 2 + 2a cos uj — 2 sin 2 uo 



(1) 



where a is the parameter of the allpass filter A(z): 



A(z) = 



a + z 



(2) 



1 + az- 1 ' 

Figure || shows the mapping functions calculated for 
some negative values of the parameter of the allpass. 
The more negative is a, the more warped are the 
modes, especially in the low frequency range. This 
result has an intuitive interpretation if we consider the 
phase delay of the allpass: the more negative is a, the 
more delayed by the allpass filter are the lower fre- 
quencies traveling into the mesh respect to the higher 
ones. 

Choosing a = —0.45 (curve in the middle), the warping 
so introduced limits the modal dispersion to very low 
values. Figure |^ shows the frequency response of the 
same TWM simulated before, after warping. Using the 



Figure 4: Impulse response taken at the center of a 
warped TWM (size 24 x 24) excited by an impulse 
at the same point. Theoretical positions of the odd 
modes resonating in a mebrane, rescaled to align the 
fundamentals (o). Positions of the same modes affected 
by residual dispersion (x). 

same notation of figure |2|, the response is compared 
with the ideal positions of the modes, rescaled to align 
the fundamentals (o), and with the real positions of 
the same modes, affected by the residual dispersion 
( x ) . The improvement in terms of precision in the 
alignment of the modes is evident by comparison of 
crosses and circles in figures § and [|. 

3 COMPUTATIONAL PERFORMANCE 

Figure || shows a plot of the dispersion factor after 
warping. Dispersion is below 2% in a range equal 
to 75% of the whole band, then it climbs to the 
maximum. From a perceptual viewpoint, it is not 
clear how much tolerance we might admit in the fre- 
quency positions of high order partials of a drum. Fre- 
quency deviation thresholds should be derived from 
subjecti ve experimentation, as it wa s done for piano 
sounds | Rocchesso and Scalcon, 1999 1. Certainly, fig. § 
shows an unnatural compression of modes that results 
in a decrease in brightness. Moreover, the frequency 
distribution of resonances brings us some information 
about the shape of the resonating object, for instance 
making it possible to discriminate a circular from a 
square drum. The warped TWM preserves the cor- 
rect distribution of modes quite well up to 75% of the 
frequency band. 

By comparison of figures [l] and we can note that a 
similar precision is achieved by a TWM if the waveg- 
uides are reduced to one third of the original length 
(the mesh is nine times denser). Then, the fundamen- 
tals can be aligned in the two models by multiplying 
times 1.75 the sampling rate of the warped TWM. Fi- 
nally, both the output signals must be lowpass filtered. 
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Figure 5: Plot of the dispersion error versus temporal 
frequency magnitude in the warped TWM. 
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Mult 


Memory 


TWM 


99 


9 


54 


WTWM 


40.25 


22.75 


22.75 


FDS 


54 


9 


18 


WFDS 


17.5 


8.75 


7 



Table 1: Performance of the TWM (FDS) vs. warped 
version ('W) in terms of sums, multiplies and mem- 
ory locations. Both the TWM (FDS) and its warped 
version allow the same dispersion tolerance. 

From these considerations, a comparison of the TWM 
versus its warped version in terms of needed sums, 
multiplies and memory locations, based on dispersion 
tolerance, can be summarized in Table |l|, where the 
warped models are labeled with the prefix W', and the 
allpass filter is supposed to be implemented in canon- 
ical (2 multiplies, 1 delay) form. Both the straight 
mesh and FDS implementations are considered. The 
number of multiplies can be further reduced at the ex- 
pense of more memory by using one-multiply allpass 
filter structures. 

4 CONCLUSION 

A new technique to reduce modal dispersion in a wide 
frequency range in TWM and triangular FDS models 
of 2D resonators has been presented. This technique 
is based on first-order allpass filters embedded in the 
mesh, and it requires an increase in temporal sampling 
rate accompanied by lowpass filtering of the output 
signal. The resulting warped TWM is shown to be less 
expensive in terms of computing resources and memory 
consumption than oversizing a TWM or FDS model. 
The coefficient of the embedded allpass filters is also a 



parameter that can be controlled to introduce tension 
modulation or other more exotic effects. 

References 

[Aird et al., 2000] M. Aird, J. Laird, and J. Fitch. 
A Model for a Drum Using Interpola- 
tion/Deinterpolation Techniques to Interface 
2-D and 3-D Waveguide Meshes In Proc. Inter- 
national Computer Music Conference, elsewhere 
in these Proceedings, Berlin, Germany, 2000. 

[Fontana and Rocchesso, 2000] F. Fontana and 
D. Rocchesso. Signal-theoretic characterization 
of waveguide mesh geometries for models of two- 
dimensional wave propagation in elastic media. 
IEEE Trans. Speech and Audio Processing, 2000. 
In press. 

[Fontana and Rocchesso, 1998] Federico Fontana and 
Davide Rocchesso. Physical modeling of mem- 
branes for percussion instruments. Acustica, 
83(l):529-542, 1998. 

[Harma, 1998] A. Harma. Implementation of recursive 
filters having delay free loops. Proc. IEEE Int. 
Conf. Acoustics, Speech, and Signal Processing, 
vol. 3, pp. 1261-1264, Seattle, Washington, 1998. 

[Harma, et al., 2000] A. Harma, M. Karjalainen, L. 
Savioja, V. Valimaki, U. K. Laine, and J. 
Huopaniemi. Frequency-warped signal process- 
ing for audio applications. AES 108th Int. Conv., 
preprint no. 5171, Paris, France, 2000. 

[Moorer, 1983] James A. Moorer. The Manifold Joys 
of Conformal Mapping: Applications to Digital 
Filtering in the Studio. Journal of the Audio Eng. 
Soc, 31(ll):826-840, 1983. 

[Savioja, 2000] L. Savioja and V. Valimaki. Reduc- 
ing the dispersion error in the digital waveguide 
mesh using interpolation and frequency-warping 
techniques. IEEE Trans. Speech and Audio Pro- 
cessing, 8(2):184-194, 2000. 

[Rocchesso and Scalcon, 1999] Davide Rocchesso and 
Francesco Scalcon. Bandwidth of perceived in- 
harmonicity for physical modeling of dispersive 
strings. IEEE Trans. Speech and Audio Process- 
ing, 7(5):597-601, 1999. 

[Strikwerda, 1989] J. Strikwerda. Finite Differ- 
ence Schemes and Partial Differential Equations. 
Wadsworth & Brooks, Pacific Grove, CA, 1989. 

[Van Duyne and Smith, 1993] Scott A. Van Duyne 
and Julius O. Smith. The 2-D digital waveguide 
mesh. In Proc. IEEE Workshop on Applications 
of Signal Processing to Audio and Acoustics, Mo- 
honk, NY, 1993. 



